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ABSTRACT 

Studying the radial variation of the stellar mass function in globular clusters (GCs) has proved a valuable tool 
to explore the collisional dynamics leading to mass segregation and core collapse. Recently, Pasquato et al. 
(2009) used the mass segregation profile to investigate the presence of an intermediate-mass black hole (IMBH) 
in NGC 2298. As a relaxed cluster with a large core, M 10 (NGC 6254) is suitable for a similar investigation. 
In order to study the radial dependence of the luminosity and mass function of M 10, we used deep high 
resolution archival images obtained with the Advanced Camera for Survey (ACS) on board the Hubble Space 
Telescope (HST), reaching out to approximately the cluster's half-mass radius (r/,m), combined with deep Wide 
Field and Planetary Camera 2 (WFPC2) images that extend our radial coverage to more than 2rhm- From our 
photometry, we derived a radial mass segregation profile and a global mass function that we compared with 
those of simulated clusters containing different energy sources (namely hard binaries and/or an IMBH) able 
to halt core collapse and to quench mass segregation. A set of direct N-body simulations of GCs, with and 
without an IMBH of mass 1 % of the total cluster mass, comprising different initial mass functions (IMFs) and 
primordial binary fractions, was used to predict the observed mass segregation profile and mass function. 

The mass segregation profile of M 10 is not compatible with cluster models without either an IMBH or 
primordial binaries, as a source of energy appears to be moderately quenching mass segregation in the cluster 
Unfortunately, the present observational uncertainty on the binary fraction in M 10 does not allow us to confirm 
the presence of an IMBH in the cluster, since an IMBH, a dynamically non-negligible binary fraction (^ 5%), 
or both can equally well explain the radial dependence of the cluster mass function. 

Subject headings: globular clusters: general — globular clusters: individual(M 10) — stars: luminosity func- 
tion, mass function 



L INTRODUCTION 

Globular Clusters (GCs) are important astrophysical labo- 
ratories for the study of both stellar evolution and stellar dy- 
namics. In recent years it has become clear that these two as- 
trophysical processes cannot be studied independently: phys- 
ical interactions between single stars as well as the formation, 
evolution, survival and interactions of binary systems have a 
significant role i n the evolution of the clusters and of their stel- 
lar populations (ISpitzeiill987h . Such interactions change the 
energy budget of the cluster and therefore influence the time 
scales on which mass segregation, core collapse and other dy- 
namical processes occur. On the evolutionary side, they can 
generate objects that cannot be explained by standard stel- 
lar evolution (lik e blue st ragglers. X-ray binaries, millisecond 
pulsars, etc.; see lFerrarol l2006 . and references therein). 

Recently, various clues have emerged that point to the 
presence of intermediate-mass black holes (IMBHs) in 
GCs. On the theoretical side at least three different for- 
mation mechanisms have been proposed, namely merg- 
ing of stellar-mass black holes via four-body interactions 
dMiller & Hamilton I2002h. runaway merging of massive stars 
(|Portegies Zwart et a 1. 1 120041) . and the death of Pop III stars 
in the early universe, which would leave IMBHs behind as 
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remnants, as shown e.g. by Madau & Rees I (|2001") (see 
also Trenti & Stiavelli 2007a). Detailed collisional N-body 
simulations (Baumgardt et al. 2005; Trenti et al. 2007b) and 
theoretical arguments (Heggie et aL ,200 7) have shown that 
an IMBH can be at the origin of a shallow central cusp 
in the surface brightness profile (SBP) of GCs that have 
bee n observed in the projected density profile of a f e w GC s 
by Novola & GebhardtJ (120061 120071) . iLanzoni et alJ (HqSI, 
and Ibata et al. (2002). 

Line-of-sight v elocity st udies were undert ak en b y 
iBaumgardt et al.' ('2003a^, Ivan den Bosch et al. I (120061). 
and Chakrabartv (2006) on M 1 5, by 
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(120021) . iBaumgardt et al. (2003b), Gebhardt et al. (20051 
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on G 1, by^ oyola et al. (2008), Anderson & van der Marel 
( 120091) . ISoUima et all (120091) on uj Cen, resulting f o r no\y 
in no undisputed definitive detection. Ilbata et al] ( l2009h 
recently found a stellar density cusp and a velocity dispersion 
increase in the center of the globular cluster M 54. This could 
be explained by the presence of a ~ 94OOM0 black hole if 
the cusp stars possess moderate radial anisotropy. Proper 
motion studies from Hubble Space Telescope (HST) data 
are likely the best way to attack the problem, but they are 
technologically challenging and inherently time-consuming, 

requiring multi-epoch HST-qua lity observations . 

Numerical simulations by Baumgard t et al. I ([2004) have 
shown that IMBHs produce a clear photometric signature in 
GCs, in terms of a larg e core to half-mass-radius ratio. More- 
over, iGill et alH ( 12008) show that a quenching of mass seg- 
regation is predicted by N-body simulations of clusters har- 
boring an IMBH, resulting in an observable signature on the 
radial dependence of the stellar mass function in GCs, which 
can be used to either suggest or rule out the presence of an 
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IMBH. IPasquato et all (l2009h provide the first application of 
this method, using the observed mass segregation profile to 
argue against the presence of an IMBH in NGC2298. More 
generally, studying the radial dependence of the stellar mass 
function in collisionally relaxed GCs can lead to valuable in- 
sights on the underlying dynamics of the relaxation processes, 
such as stellar evaporation, mass segregation and core col- 
lapse. 

Mass segregation occurs in star clusters as two-body re- 
laxation processes drive the system towards energy equipar- 
tition. Heavier stars sink towards the center of the cluster, 
while the lighter stars preferentially live in the outer regions 
and are more likely to attain escape velocity from the system. 
If an additional energy source is present, such as a dynami- 
cally significant population of binaries or an IMBH, the mass 
segregation and the core collapse processes are slowed down 
and the arising gradient of the IMF is less steep. Also a sig- 
nificant population of stellar mass black holes can in princi- 
ple segreg ate in the center of the cluster and sustain a large 
core (see iMackev et al.l 120081) . thereby acting as an energy 
source . This effect is observed in our N-BODY simulations 
with a lSalpeted (llQSS *) IMF (see Section ID, which produces 
comparatively more remnants than a Miller & Scalo (1979) 
IMF. However, in such simulations this effect typically lasts 
a few Gyr, because stellar mass black holes kick each other 
out of the system via dynamical interactions, and is there- 
fore unlikely to be at work in M 10, which is approximately 
11.8 ± 1 . 1 Gyr old (Salaris & Weiss 2002). In any case, a 
ISalpeted (119551) IMF likely overestimate the number of large 
stellar-mass black holes produced in a real GC, so any effect 
of such a dark core on the long-term dynamics of the cluster 
is effectively ruled out. 

IMBHs and binaries both work towards reducing the degree 
of mass segregation attained by a GC at a given (dynamical) 
age. A realistic binary fraction can partly mimick the effects 
of an IMBH, so a detailed knowledge of the binary fraction 
allows more robust conclusions either in favor or against the 
presence of an IMBH. So far, however, only for a few GCs 
has a reliable binary fraction been estimated with robust pho- 
tometric methods based on main sequence ( MS) broadening 
(see ISolhma et al.ll2007t IMilone etaH 12008. and references 
therein), so the type of study that we discuss here is only pos- 
sible for a handful of clusters. 

M 10 is a reasonable candidate for harboring 
an IMBH since it displays a ratio of almost 0.4 
(iMcLaughlin & van der Marell 120051) between the core 
radius and the half-light radius, which in dynamically old 
clusters is usually interpreted as the si gnature of an extra 
energy source (see IVesperini & Chernoffi,19 94). In addition, 
MIO is an ideal candidate for the mass segregation method 
described and employed for NGC2298 by Pasquato et 
al. (2009) because , with a mass of about 1.5 x lO'' M© 
dMcLaughl in & van d er Marel 2005J), it is dynamically re- 
laxed (log(f/j) = 8.863^5; Harris 1996) and has not been overly 
influenced by the Galactic tidal field (r,/r/,„, = 11.86 and Rgc 
= 4.6 kpc; Harris 1996). These conditions ensure that it has 
experienced many initial half-mass relaxation times and thus 
has achieved equilibrium with respect to mass segregation. 

The structure of the paper is as follows. We present the 
data analysis in Sect. [2] a nd the resulting color-magnitude di- 
agram of M 10 Sect. 12.21 The radial dependence of the mass 
and luminosity function is showed in Sect. |3] while in Sect. |4] 
the global mass function (GMF) of the cluster is derived. In 



TABLE 1 
Observations 



Camera 


Filter 


#exp. 


Exp. Time 


ACS 


F606W 


4 


90s 




F814W 


4 


90s 


WFPC2 


F606W 


1 


1100s 




F606W 


3 


1200s 




F814W 


1 


1100s 




F814W 


3 


1200s 



Sect. |5] we compare the observed mass segregation profile 
with predictions from our N-body simulations. A summary 
and conclusions follow in Sect.|6] 

2. OBSERVATIONS AND ANALYSIS 

The photometric data used here consist of two main data 
sets. The "ACS sample" consists of a set of HST images 
of the core region obtained with the Advanced Camera for 
Survey (ACS) through the F606W and F816W filters (GO- 
10775; P.I.: A. Sarajedini; see Table©. The ACS Wide Field 
Channel (WFC), employs a mosaic of two 4096 x 2048 pixels 
CCD, providing a plate-scale of 0"05/pixel and a total field 
of view (FOV) of 3f4 x 3. '4. The dithering pattern used in 
the observations allows us to sample cluster stars from the 
very center out t o the nominal half-mass radius (r/, ^ 1.'81; 
[Harris et"an (11996.) ). Since we were interested in deriving ac- 
curate photometry of the cluster's low-mass MS population 
(M < O.BMq), we retrieved from the ESO/ST-ECF Science 
Archives only images with long exposure times (fexp = 90 s). 
All the flat-fielded images (FLT) were properly corrected for 
geometric distortions and for the effective a rea of each pixel, 
following the prescriptions of lSirianni et alj (12005)- 

The "WFPC2 sample" consists of a series of deep HST im- 
ages obtained with the F606W and F814W filters of the Wide 
Field Planetary Camera 2 (WFPC2) on 1995, September 30 
during Cycle 5 (Prop.: 6113; PL: Paresce). A description of 
the number of exposures and their duration is given in Table[T] 
The WFPC2 field is located about 3' South- West (SW) of the 
cluster center Also in this case the data were obtained through 
the ESO/ST-ECF Science Archives and were pre-processed 
through the standard HST pipeline. We did not reduce the PC 
chip, since its small FOV (^ 40" x 40"), combined with the 
extremely low stellar density at this distance from the cluster 
center (~ 4'), sampled a very small number of objects, with 
no appreciable improvement to the general statistic. A map of 
the entire data-set and is shown in Figure[T] 

2.1. Data reduction 

The goal of this work is to obtain accurate photometry of 
MS stars down to the lowest mass allowed by the exposures. 
Since M 10 is a reasonably sparse GC (c = 1 .4C0; Harris et al] 
(119961) ) the images are not affected by severe crowding, even 
in the very central regions. The photometry of two data-sets 
was performed as follows. We used DAOPHOTII/ALLSTAR 
(IStetsonlll98 7) to find stars, derive a spatially varying point 
spread function (PSF) and perform the first measurements for 
all stars in every single image. We then used DAOMATCH 
and DAOMASTER to match all stars in each chip, regardless 
of the filter, in order to find an accurate coordinate transforma- 
tion between the frames. The matched solutions were then fed 
to M0NTAGE2 in order to build a stacked image of each chip. 

* Central concentration, c = log{rt /vc) 
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In this way we could eliminate all the cosmic rays and ob- 
tain an image of the stars with the highest signal-to-noise ra- 
tio. We ran the DAOPHOT/FIND routine and the PSF-fitting 
on the stacked image in order to obtain a deep master star 
list. The mast er list was then used as input for ALLFRAME 
dS tetsonl l 19941 ). which simultaneously determines the bright- 
ness for stars in all frames while enforcing one set of cen- 
troids and one transformation between all images. Finally, all 
the magnitudes for each star were normalized to a reference 
frame and averaged together, and the photometric error was 
derived as the standard deviation of the repeated measures. 

We note here that the presence of saturated stars in these 
deep images is the major source of artefacts causing spuri- 
ous detections, especially in the core regions where most of 
the red giant stars are located. In order to clean our cata- 
logues from false detections, we used a statistical approach 
similar to the one proposed by ICool et al.l (11996 ). that em- 
ploys the sharpness (sh) parameter provided as an output by 
ALLFRAME as the source quality diagnostic. By plotting 
sh as a function of the magnitude we found that spurious ob- 
jects, detected around the haloes and along diffraction spikes 
of saturated stars, have sharpness values that differ from the 
ones representative of bona-fide stars, which typically have 
-0.15 < sh < 0. 15. In Figure|2]we show that the effectiveness 
goodness of the selection in sharpness (left panel) in identi- 
fying false detections is confirmed by the fact that the objects 
with large sharpness values have also larger photometric er- 
rors (right panel). 

We finally transformed the instrumental F606W (V) and 
F814W (I) magnitu des to the VE GAMAG s ystem following 
the pre scriptions of ISirianni et al.l (|2P05i) and iHoltzman et alj 
(Il995lfl for the ACS and WFPC2 samples, respectively The 
relative star coordinates of both samples were transformed to 
the absolute right ascension and declination values (J2000) 
using the wide field catalogue published by Polla rd et aL 
( 120051) as secondary astrometric standard catalogue. 

2.2. The color magnitude diagram 

In Figure [3] we show the color-magnitude diagram (CMD) 
of the entire data set. The WFPC2 final catalogue (left 
panel) contains 4995 objects. The same data se t have al- 
ready been reduced bylDe Marchi & Paresc3 (11996 ) and later 
bv .Piotto & Zoccalil Il99'9l) . By comparing our photometry 
with the CMD of the latter (see their Figure 1), no differences 
are evident. The consistency of the two photometric analyses 
is further confirmed by a direct comparison of their published 
luminosity function (LF; see Table 4) and the one calculated 
in this work for the WF2 chip (see Section [3] for a detailed 
discussion). 

Due to the long exposure time of the WFPC2 images, the 
stars at the turn-off (TO) level in this data set are saturated. 
In order to obtain a CMD in the WFPC2 area with a I mag- 
nitude range comparable with that of the ACS data set, we 
decided to complement our observation s with a V vs. iV -I) 
ground based catalogue from Ro senbergetaH (l2000fl. This 
data set was obtained using the 1.0 m Jacobus Kapteyn Tele- 
scope (JKT) at La Palma (Canary Island). We transformed the 
standard Johnson V, / magnitudes of the JKT catalogue to the 
WFPC2 VEGAMAG system adopting color term as derived 

^ Notice that we refeiTed to Section 5.1 of WFPC2 Data Handbook for 
updated values of Photometiic Zeropoints 

* The photometric catalogue can be retrieved at the web page 
[http://www.astro.unipd.it/globulars/ 



by comparison of all the non-saturated stars in the WFPC2 
chips that were also measured in the ground-based images. 
No color term is necessary for the / band. The uncertainties 
in the V zero points and in the (V-I) color transformations are 
of the order of 0.007 mag. Due to the extremely low crowd- 
ing conditions of the area covered by the WFPC2 and to the 
relative proximity of M 10 {R^^n = 4.4 kpc; Hams 1996), the 
completeness of the ground based catalogue in this area is ex- 
pected to be very high. By comparing the observed LFs from 
the WFPC2 and JKT data in the same area, we found a good 
agreement in the magnitude range 19 < / < 20 (corresponding 
to the mass range of approximately O.6M0 < M < 0.7Mq). 
We therefore decided to use the JKT catalogue for stars with 
/ < 20 and the WFPC2 catalogue at fainter magnitudes. 

In the right panel of Figure [3] we show the CMD as derived 
from the photometric reduction of the ACS data set. The fi- 
nal catalogue contains 56812 objects. The MS of the cluster 
is clearly visible and well sampled well above the TO that 
occurs at I'^ 17.5. A population of candidate blue stragglers 
departing from the TO along a bright extension of the MS is 
also recognizable. For / < 16 saturation starts to occur, thus 
making the photometry of brighter stars unreliable. 

According to Piotto & Zoccali (1999), due to the Galac- 
tic latitude of MIO {b = 23°), the contamination from fore- 
ground/backgrou nd stars is small and should not effect the LF. 
On the other hand lPoUard et al.l (l2005h . using wide field imag- 
ing of evolved stellar populations, claim that field contami- 
nation, mainly from the Galactic disk, should be taken into 
account. Since no direct measurement of the field contami- 
nation in our magnitude range is possible, we used a statisti- 
cal ap proach. A c atalo gue of stars using the Galactic model 
from Robin et aTl (120031) coveringan area of 1 square degree 
around the cluster center was usec0. After scaling for the ACS 
and WFPC2 FOV, we found a contamination of 425 and 159 
stars in our magnitude range respectively for each sample. 
This means that, even in the worst case of the WFPC2 sam- 
ple, where the number of sources is small, field contamination 
would be of the order of ^ 3 %, and thus can be neglected. 

The MS mean ridge lines for the two data sets (open circles 
in Figure O were computed by using a 2nd order polynomial 
to fit the locus of MS stars, after rejecting those farther than 
2 a from the best fit line, where a is the combined photometric 
uncertainty in V and /. Following De Marchi & Pulone (2007; 
hereafter DP07), we applied this same d-clipping approach 
to identify the bona-fide stars to be used in computing the LF 
and MF of the cluster and, in order to increase the statistics, 
we decided to accept as bona-fide stars all those within 2.5 cr 
of the ridge line of best fit. The catalogue obtained after this 
procedure, containing 46 407 and 4 390 stars respectively in 
the ACS and WFPC2 field, will be used hereafter to derive 
the LF and MF. We mark as shaded regions in the CMDs of 
Figure [3] the limits below which the photometric complete- 
ness, as detailed in Section 2.3, falls below 50 %. All the stars 
in these region are not used for the LF and MF calculation. 

The MS ridge line that we derived in this way agrees very 
well with the models of Baraffe et al. (1997). The level of 
agreement between the models and the observed data can be 
fully appreciated in Figure [3] where the models for metallic- 
ity [M/H] =-1.3 (dashed line) and [M/H] = -1 .0 (soHd line) 
are sh own. As noted by DP07, the models of Baraff e et alJ 
( Il997h are calculated for the F606W and F814W filters of the 
WFPC2, which are slightly different from those on baord the 

' http://model.obs-besancon.fr/ 
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ACS. We used the same method proposed in their work in 
order to translate the fiher from the WPFC2 to the ACS sys- 
tem. Briefly, by using synthetic model atmosphere from e.g. 
the ATLAS9 library of Kurucz (1993), it is possible to cal- 
culate the magnitude difference in the same filter for the two 
cameras as a function of the effective temperature of the star. 
As shown in Figure [3] the agreement of the translated model 
in the ACS sample (right panel) is very good. The best fit 
of the model to the fiducial line corresponds to metallicity 
[M/H] =-1.0, distance modulus (m-M) = 14.21 and color 
excess E{B-V) = 0.26. These values are in excellent agree- 
ment with those of lPiotto & Zoccalil ( 119991) . namely (m-M) = 
14.20 and E(B-V) = 0.29. The metallicity of the best fitting 
model, [M/H] = -1.0, is fully consistent with the va lue mea- 
sured observationally by ICarretta & GrattonI (Il997h . namely 
[Fe/H] = -1.41, whereas those for [M/H] = -1.3 deviate 
in colo r in the lowe r part o f the CMD, as already pointed 
out by iPoUard et al.l ( |2005|) . We used an isocrone calcu- 
lated for_an_ageof lOGyr, the only available at the moment 
from Baraffe et alj ([1997,). It is clear from their Figure 3a that 
the age plays a minor when void effect in the range of masses 
considered in this work. In light of the good agreement be- 
tween the observed ridge line and the models of Baraffe et d] 
(11997), we feel confident that we can use the latter to translate 
magnitudes into masses. 

2.3. Photometric completeness 

A crucial step in determining an accurate LF is the estima- 
tion of the photometric completeness of the data. Complete- 
ness corrections were determine d via artificial s t ar test s fol- 
lowing the method discussed bv lBellazzini et al] (120021) . We 
first generated a catalogue of simulated stars with a V mag- 
nitude randomly extracted from a LF modeled to reproduce 
the observed LF in the V band. Then the / magnitude was 
assigned at each sampled V magnitude by interpolating the 
mean ridge line of the cluster (see Figure[3]l. The coordinates 
of the simulated stars where calculated in the coordinate sys- 
tem of the reference frame and then translated to each single 
image using the same matching solutions found in the data 
reduction phase (see Section lZTT i. 

It is crucial to avoid interference between the artificial stars, 
since in that case the output of the experiments would be bi- 
ased by artificial crowding not present in the original frame. 
For this reason the frames were divided in a grid of cells of 
fixed width (^ 15 pixel, i.e. more than 5 times larger that 
the mean FWHM of the stars in the frames) and only one star 
was randomly placed in such a box in each artificial test run. 
The artificial stars were added to the real images using the 
DAOPHOT/ADDSTAR routine. The reduction process was 
repeated on the artificial images in exactly the same way as 
for the scientific one s and applying the same selection criteria 
described in Section lTTI More than 100000 stars were even- 
tually simulated in each chip, for a total of more that 500,000 
stars in the entire data set. The photometric completeness (/) 
was then calculated as the ratio of the number of stars recov- 
ered after the photometric reduction {Nout) and the number of 
simulated stars (A^,„). / is expressed as a function of / magni- 
tude. 

Due to the strong gradient in stellar density from the cen- 
ter of the cluster outwards, there is a significant dependence 
of the photometric completeness with radial distance, even 
though M 10 is a reasonably sparse GC and the images do 
not suffer from dramatic crowding. Nevertheless, in order 
to take this radial effect into account, we divided the en- 



tire FoV in 6 zones characterized by a similar completeness. 
In Figure [T] the black full circle show the regions numbered 
from 1 to 6. The innermost four regions are included in 
the ACS FOV and are concentric annuli all centered on the 
nominal cluster center. We notice here that the center used 
in this work has been determined from accurate star counts 
of resolved stars in the central regions of the cluster from 
archival WFPC2 imaging of M 10 (Prop. 6607, PI. F Fer- 
raro), as explained in a forthcoming paper (Dalessandro et 
al. 2010; in preparation). The adopted center of gravity is 
aj20oo = 16*'57'" 8^92 , 5nm) = -4° 05' 58"07, which is in full 
agreement with the one of Harris et al. ( 1996). 

Region 5 includes both the WF3 and WF4 chips of the 
WFPC2, which cover and area of equal crowding condition 
at the same radial distance, whereas Region 6 corresponds the 
whole WF2 chip. In Figure |4] we report the resulting photo- 
metric completeness of each of the six regions as a function 
of the /-band magnitude. It is important to notice that, even in 
the very central region of the cluster where saturation biases 
the completeness, artificial star tests show that we sample the 
MS stellar population of the cluster down to ^ 5 mag below 
the TO (see solid line in lower panel), with a completeness in 
excess of 50 %. 

3. THE LUMINOSITY AND MASS FUNCTION 

In order to derive the LF, we selected from our bona-fide 
photometric catalogies all stars within each of the six regions 
mentioned above and sorted them as a function of the /-band 
magnitude. The resulting observed LFs (OLFs) are shown as 
full circles in Figure|5]and as filled circles after correction for 
photometric completeness. The observed and completeness- 
corrected LFs are also provided in Table|2] In Figure[5]we also 
show, as a dashed line, the LF published by Piotto & Zoccall 
(1999) in the WF2 frame, which appears in excellent agree- 
ment with ours. 



The solid line shown in Figure |5] are the theoretical LFs 
(TLFs) obtained by multiplying a simple power-law mass 
function (MF) of the type dN /dm oc by the derivative 
of the mass-luminosity relationship of Baraffe et al. (1997) 
for [M/H] = -1.0. We have adopted a distance modulus 
{m-M) = 14 .21 and a color excess E{B-V) = 0.26 as found in 
Section l2!2l The power-law indices that correspond to the best 
fitting models for the six regions, from Region 1 to Region 
6, are respectively a = 0.7,0.4,0.1,-0.3,-0.6,-0.9. With the 
notation used there, the Salpeter IMF would have a = -2.35 
and a positive index implies that the number of stars is de- 
creasing with decreasing mass. The wide range of MF slope 
indices from the cluster core to 2.5 rhm indicates that stars in 
M 10 have experienced the effects of mass segregation. Mas- 
sive stars are much more segregated in the central regions 
were the index is positive. The MFs are almost flat in the 
intermediate regions with an inversion of the trend outside, 
where the number of low mass stars is increasing compared 
to the more massive ones. We address this issue more quanti- 
tatively in Section|4] 

4. THE GLOBAL MASS FUNCTION 

In order to further investigate the dyn amical state of 
M 10, we used the approach developed by iGunn & GriffinI 
(T97§). We ran the multi-mass Michie-King (MK) code 
originally developed by .Mevlaa (.1987, .1988.) and modified 
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TABLE 2 

Luminosity functions measured in each of the six regions. For each region, the table gives, as a function of the m^i^ magnitude, 

THE NUMBER OF STARS PER HALF-MAGNITUDE BIN BEFORE (No) AND AFTER (N) COMPLETENESS CORRECTION AND THE UNCERTAINTY ((Tn) ON N. 







Area 1 

N 


(Tn 


A^o 


Area 2 

N 


(Tn 




Area 3 
N 


(Tn 




Area 4 
N 


(Tn 


A^o 


Area 5 
N 


(Tn 


A'o 


Area 6 
N 


ctn 


17.25 


121 


131 


24 


259 


271 


30 


269 


279 


27 


511 


517 


32 


28 


28 


7 


10 


10 


4 


17.75 


195 


207 


31 


414 


434 


37 


497 


511 


37 


948 


961 


45 


55 


55 


10 


10 


10 


4 


18.25 


253 


285 
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by II\ilone et alj ( 119991) . to assess whether the observed MF 
variations with radius are consistent with the degree of mass 
segregation expected from energy equipartition due to two- 
body relaxation. In short, every model run generates a power- 
law MF of variable index a characterized by four cluster 
structural parameters i.e. the core radius (rc), the scale veloc- 
ity (Vj), the central value of the dimensionless gravitational 
potential Wo and the anisotropy radius (ra). The first selection 
of the models is done by excluding all those regions in the pa- 
rameter space that do not fit simultaneously the observed SBP 
and the central value of the velocity dispersion (cTp). 

The SBP used here was derived from accurate star counts 
obtained by combining a photometric catalogue sampling 
cluster stars in the central regions from WFPC2 high resoution 
images and multi-band wide field images from the Wide Field 
Imager (WFI) mounted on the MPI 2.2 m telescope at La 
Silla (Chile), extending out to the tidal radius {r, = 21f48; 
Har ris et al. 1996). Details about the SBP can be found in 
iDalessandro et alj (I2009t and references therein). As regards 
the velocity dispersion, the value of a,, = 6.60 Km/s comes 
from lMcLaughlin & van der Marell (l2005h . 

While forcing the MK model to reproduce the SBP and the 
observed value of the velocity dispersion will set some con- 
straints on the cluster's structural parameters and on the range 
of indices of the MF, in order to find a unique solution the 
additional condition that the model LFs agree simultaneously 
with the OLFs at different radii must be imposed (see e.g. De 
Marchi et al. 2006; DP07). This sets very stringent constraints 
on the shape of the present global MF (GMF), i.e. the MF of 
the cluster as a whole. 

We do so in Figure|6] where we compare the LFs (corrected 
for completeness, solid circles) with the theoretical MFs pre- 
dicted by the MK model at the specific radial locations, folded 
through the M-L relationship, exactly as we did for the model 
power-law MFs in Figure |5] Here, however, the shape of the 
local MF is not free, but is completely set by the two-body 
relaxation process acting on the GMF, as per the MK model. 
The best fitting model also reproduces very well the SBP with 
a core radius r^. = 50"4 and a dimensionless central potential 
Wo = 7, in excellent agreement with the values of r^. = 49 "41 



and Wo = 6.50 from lMcLaughUn & van der Marell (l2005h . 

The GMF index of the best fitting model is a = -0.7 (as 
already mentioned in our formalism a Salpeter IMF would 
have a = -2.35) in the range 0.25- 0.8 Mq. In the case of 
NGC2298 (DP07) and NGC6218 (De Marchi et al. 2006), 
an excellent match was found between the index of the GMF 
and that of the local MF as measured at the half mass radius. 
It is expected that the actual MF near the half mass radius 
should be only marginally affected by mass segregation and 
should thus be representative of the GMF of the cluster (see 
De Marchi et al. 2000, and references therein). It is interest- 
ing to note that also in the case of M 10 the index of the GMF 
(a = -0.7) is in good agreement with the value a ^ -0.6 (see 
Section O measured near the half mass radius (r/,„, = 124"; 
see Section|5]l. 

The good match between the observed LFs and those from 
the MK model suggests that, in its present state, M 10 is con- 
sistent with a condition of energy equipartition. On the other 
hand, since the MK approach provides by design a snapshot 
of the cluster evolution, it is not possible to infer from it the 
past dynamical history of the cluster, and in particular whether 
internal energy sources such as binary stars or an IMBH are 
playing any role in the evolution of M 10. However, this is 
possible using dynamical N-body models that incorporate sin- 
gle and binary star evolution; these will be discussed in the 
next section. 

5. N-BODY SIMULATIONS 

We use state-of-the-art direct N-body simulations t o model 
the dy namics of M 10. They are based on NBODY6 ( AarsethI 
l2003h . modified as discussed in iTrenti et al. I (l2007bl) to im- 
prove accuracy in the presence of an IMBH. Interactions be- 
tween particles are treated using regularization of close grav- 
itational encounters with no softening. This choice is mo- 
tivated by the need to realistically model the collisional dy- 
namics within the sphere of influence of the IMBH. An in- 
depth discussion of th e numerica l simu lations and the soft- 
ware used is found in iGill et alH ( 120081) . We use the same 
set of runs presented in that paper, with the addition of fur- 
ther simulations containing primordial binary fractions equal 
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to 1%, 3% and 5%. These primordial binaries were assigned 
binding energies in the range from e,„„, to 133e„„„, where 
emm = ('w(0))cre(0)^, adO) is the initial velocity dispersion, and 
(m(0)) the initial average mass of the svsterrPl This choice 
makes our binaries hard, with separations typically lower than 
10 AU. Binaries are likely to have a non-negligible impact on 
the results we obtain, so the addition of these runs is an im- 
portant improvement over lPasquato et akl (120091) . The masses 
of stars in binary systems have been selected independently 
according to the chosen IMF after the instantaneous step of 
stellar evolution has been applied to evolve the turn-off mass 
to 0.8 Mq. We are essentially adopt ing a random pairing 
function (see lKouwenhoven et al.lf2009l) . Binary exchange in- 
teractions are frequent, especially in the core of the system, 
so by f > 5fr/!(0) (when we consider the simulation snapshots 
for the analysis presented here) more than 50% of the bina- 
ries no longer have their initial companion. Hence the de- 
tails of the pairing function are not expected to significantly 
affect the measured amount of mass segregatio n. The ini- 
tial binary eccentricity distribution is thermal (see Trenti et al.' 
I2OO8). For a typical velocity dispersion of 10 km/s, the binary 
binding-energy distri bution corresponds to a semiaxis distri- 
bution < lOAU (see iTrenti et all 120081 l20Tol) . Binaries with 
larger separ ations w'ould b e quickly destroyed by three-body 
encounters (iHeggidl 1975b while binaries with larger binding 
energy (and smaller semiaxis) than those considered here are 
generally possible only if both components are compact rem- 
nants (otherwise their apocenter separation could be less than 
the physical radius of the stars for eccentric orbits). In any 
case su ch tightly bou nd objects are essentially dynamically in- 
ert (see iTrenti et al.l l2008, for runs that consider an extended 
binding energy range). Observations of nearby star forming 
regions and the Galactic field, find a pri mordial binary frac- 
tion which is usually quite hig h (Duque nnoy & M ayor 1992; 
iKouwenhoven et al.ll2005L 12007 ) but this refers to all binaries, 
independent of their binding energy. In our case only hard 
binaries are included in the primordial binary population, as 
softer binaries would be quickly destroyed by three-body in- 
teractions, having no long-term dynamical effect. This makes 
our choice of exploring the 0%- 10% range in primordial bi- 
nary fraction justified. 

Here we summarize the initial conditions to the various 
runs we use to explore a variety of different primordial con- 
figurations, changing the IMF to check that our conclusions 
are independent of the initial configuration, and the primor- 
dial binary fraction to study its influence on our conclusions . 
The initial distribution function is a single-mass Kh^(fT966h 
model, even though a fully realistic mass spectrum is then 
used in the N-body calculations. The number of particles is 
varied from = 1 .6 • lO"* to = 3.2 ■ 10"*. If we take the aver- 
age mass of a GC star to be about O.5M0, this numbers trans- 
late to roug hly 8 X 10^ Mq and 1 .6 x 10"^ M,^. Conversely the 
mass of MIO is estimated by iMcLaughlin & van der Marell 
(j2005) to be w 1.55 x IQ^Mq. 

None of our simulations includes primordial mass segrega- 
tion, as all stars start out with the same distribution, regardless 
of mass. Recently, primordial mass segregation was proposed 
to explain the extreme depletion of the stellar mass func- 
tions observed in several GCs, which exceeds the predictions 

The choice of the number 133 follows Heg gie et all <20O6l) . where bi- 
nary binding energies ranged from 3kTc to 400A:rc, with 7c defined as the 
central temperature of the system. 



of pri mordially non-segregated simulations ( IBaumgardt et al.l 
I2OO8I) . Primordial mass segregation, in combination with an 
external tidal field, can be used to justify a much higher rate of 
depletion of low-mass stars. Moreover, Allison et al. (200^ 
find that primordial mass segregation naturally arises in clus- 
ters initialized in a clumpy, subvirial (dynamically cold) state. 
In the present work we do not explicitely include mass seg- 
regation in the initial conditions. However, in order to draw 
our conclusions on the presence of an IMBH we consider the 
mass segregation profiles of relaxed configurations only, by 
looking at times t > 5f/,. So it is the equilibrium value of mass 
segregation that matters for our argument, and we expect it 
not to be influenced by its primordial value, similarly to what 
is found for r^/ r/,„„ which settles to an equilibrium value after 
several relaxation times, irrespective of its starting value (see 

Trenti et al. 2007c; Trenti et al. 201^). 

Initial stellar masses come from either a lSalpeterl (119551) or 
iMiller & Scalo (1979) IMF. We have adopted a lower limit 
to the IMF of O.2M0 and an upper limit of lOOM© (prior 
to the stellar evolution step). The choice of the lower mass 
limit is driven by the mass of the faintest stars that are observ- 
able in our data set with a satisfactory completeness (w 50%). 
Pasquato et al. (2009) computed control runs including stars 
with masses down to OAMq, obtaining the same results as 
in the case of runs with masses truncated at O.2M0. Stellar 
evolution is carried out in a single step to a turn-off mass of 
Mt.o. = O.8M0, at the beginning of the simulation. The details 
and th e motivation of this choice are explained by Gill et aO 
(2008). Half of the runs included in Pasqua to et al. I (l20oC 
contained primordial binaries, with a binary fraction of 10%. 
This includes runs containing both an IMBH and a 10% bi- 
nary fraction, which allow us to study the combined effect 
of the two ingredients. Qualitatively, the quenching effects 
of binaries and of an IMBH sum, leading to an enhanced de- 
pression of mass segregation. This fraction can reasonably be 
considered as the value represe ntative of the binary fra ction of 
a galactic globular cluster (see iBellazzini et ani2002 l). In this 
paper, new simulations without an IMBH but with primordial 
binary fractions of 1%, 3% and 5% are considered. These new 
runs have a Mill er & Scald (119791) IMF and 3.2 • 10"* particles. 
They are included to give a more fine-grained scale in binary 
fraction, to properly compare its dynamical effects with the 
observations. 

In about half of the simulations, an IMBH with mass 
MiMBH — 0.01, or ~ 1% of the entire cluster is present. Values 
of the M IMBH I Mgr. in the rang e 0.5% to 2% are expected by 
extendin g the iMagorrian et aP (i 1 996) scaling law to GCs. As 
shown by lGill et al. 1 (120081) . even rising the mass of the IMBH 
up to 3% of the cluster mass does not cha nge the amount of 
mass segregation found in the simulations. iGill et alH (l2008h 
further suggest that, in general, the dependence of the quench- 
ing of mass segregation on IMBH mass is weak, given that 
the quenching mechanism is based on the scattering of stars 
by the binary system formed in the GC center by the IMBH 
and a massive companion (e.g. a stellar-mass black hole). 
The tidal field from the parent galaxy is taken into account 
in the simulations assuming circular orbits with the tidal cut- 
off radius s elf-consistent with th e initial cluster concentration 
parameter. Dinescu et alJ (119991) find that MIO has a moder- 
ately elliptic orbit around the galactic center, with eccentricity 
e = 0.19 ±0.05. This low value of eccentricity makes the cir- 
cular orbits adopted in this work a satisfactory approximation 
to the orbit of MIO. For further details of the tidal force treat- 
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ment see lTrenti et al. I (l2007cl) . 

5.1. Comparison with the observations 

We applied the same algorithm used by iPasquato et al. I 
(|2009) on NGC 2298 to the catalogue of stars we obtained 
for M 10, in order to produce a radial mass segregation pro- 
file. A half-mass radius was obtained for main sequence stars 
(with masses ranging from 0.26 to 0.8 Mq), by computing 
their total mass within the radial range covered by the ob- 
servations and finding the radius which contains half of it in 
projection. A non-parametric mass growth-curve was con- 
structed via spline interpolation as in IPasquato et al.1 (120091) . 
We obtain r/„„ = 124", which is sli ghtly larger than the half- 
light value of 108.6" determined bv lHarris et al] (Il996h . This 
is not surprising and is per se an indication of mass segre- 
gation, since bright giants (with masses around the turn-off 
mass) dominate the luminosity profile. We used the obtained 
half-mass radius to normalize our mass-segregation indicator, 
i.e. th e mean mass of MS sta rs as a function of radius, fol- 
lowing l^sguatoeTalll (120091) . The mass segregation profile 
ATO(r) is constructed by plotting the mean mass of a main- 
sequence star (normalized by subtracting its value at the half 
light radius) as a function of radius. 

We compared the observational mass segregation profile to 
the confidence intervals predic t ed by means of the set of simu- 
lations used in lPasquato et al. I (12009 ) with the addition of the 
three new IMBH-free runs with a 1%, 3%, and 5% primordial 
binary fraction, consistently defining t he range o f observable 
MS stars from 0.26 to 0.8 Mq. iGill etlTI (|2008) show that 
the simulations reach a stationary configuration with respect 
to mass segregation after approximately 5 relaxation times. 
For the sake of comparing them to the observed mass segre- 
gation profile of M 10, the simulations were analyzed between 
7 and 9 initial-half-mass relaxation times. 

Figure |7] (left panel) shows the M 10 mass segregation pro- 
file (red line) superimposed to the green and blue shaded ar- 
eas corresponding to the Icr confidence regions defined by 
the simulations without an IMBH and with an IMBH of mass 
1% of the total cluster mass respectively. The observational 
profile lies within the superposition of the two shaded areas. 
Therefore, while it is possible that an IMBH is lurking in the 
center of M 10, we cannot rule out alternative explanations. 
In particular, the lower envelope of the green shaded area is 
defined, towards the center, by runs containing primordial bi- 
naries with an initial binary fraction of 10%. This is expected, 
because binaries act as an energy source capable of halting 
core collapse, as shown by Trenti et al. (2007b c) and con- 
sequently can partially mimick the behavior of an IMBH, so 
they give rise to a shallower profile with respect to runs con- 
taining only single stars. The right panel of Figure [T] shows 
the same plot of the left one, where only the 32k-particle sim- 
ulations with no primordial binaries were included (binaries 
do form dynamically in these runs, but only in small num- 
bers). In this case the data would be incompatible with the 
no-IMBH scenario, and we would conclude that an IMBH is 
likely present in the core of M 10. 

The two plot described above clearly show the crucial role 
played by binary stars in the dynamical evolution of the clus- 
ter. In order to investigate in more detail the effect of the 
binary fraction in the predictive capability of our model, we 
compare the observed mass segregation profile of M 10 with 
the prediction of the 32k-particle IMBH-free simulations with 
3% and 5% primordial binary fraction (respectively left and 
right panel in Figure O . It is apparent that these simulations 



can reproduce the observations fairly well, i.e. that in the ab- 
sence of an IMBH a binary fraction in this range is required 
to explain the observations of mass segregation in M 10. 

We therefore conclude that a precise measurement of the 
MS binary fraction in M 10 would allow us to reach a firmer 
conclusion on the presence of an IMBH. A dynamically sig- 
nificant fraction of binaries is a viable alternative to the IMBH 
hypothesis for explaining the observed mass segregation pro- 
file, since both can be active energy sources capable of halting 
core collapse and quenching mass segregation. 

6. SUMMARY AND CONCLUSIONS 

In this work we present a study of the dynamical state of 
the globular cluster M 10 through the characterization of the 
radial properties of its stellar LF and MF. We used a com- 
bination of archival deep ACS and WFPC2 images to sam- 
ple cluster stars from the core regions through to ^ 2.5 r/, 
in the mass range 0.25- 0.8 Mq. The full data set was di- 
vided in six regions at different distances from the cluster 
center and characterized by the same photometric complete- 
ness. The LF was calculated for each region an d converted to 
a MF using the mass-luminosity relationship of lBaraffe et alj 
(1997) for the metallicity of the cluster. Each of the lo- 
cal MFs was fitted using a simple power law of the type 
dN /dm cx m" and the best fitting indices were found to be 
a = 0.7, 0.4, 0.1, -0.3, -0.6, -0.9 from the center outwards. A 
positive value of a means that the number of stars decreases 
with decreasing mass. This radial change in the MF index is 
a clear sign of mass segregation in M 10. To better charac- 
terise this effect, we built a multi-mass Michie-King model 
to reproduce the observed radial variation of the MF and the 
cluster's SBP and veolcity dispersion and found that the dis- 
tribution of the stars in the cluster is compatible with a condi- 
tion of equipartition of energy. The GMF as obtained by the 
Michie-King model is a power-law with index a = -0.7, in 
agreement with the value a ^ -0.6 near the half mass radius 
(rh,„ = 124"). 

In order to investigate in more detail the dynamical his- 
tory of M 10, we used N-body simulations that include the 
presence and evolution of binaries, with the goal to under- 
stand whether the current mass segregation profile requires 
the presence of an IMBH. Within the current uncertainties, 
we cannot exclude an IMBH of mass ^ 1 % of the total cluster 
mass, although a plausible alternative explanation of the ob- 
served mass segregation profile is the presence of a substantial 
amount of binaries, likely primordial in origin, capable of sig- 
nificant dynamical effects. We show that N-body simulations 
initialized with a primordial binary fraction in the range of 
3 - 5 % are effectively able to correctly predict the observed 
mass segregation profile of M 10. 

MIO is a better candidate than NGC 2298 for an ob- 
servational follow-up aimed at detecting an IMBH. While 
the latter cluster does not contain an IMBH according to 
Pasquato et al. (2009), M 10 shows a shallow mass segrega- 
tion profile which could be interpreted as the fingerprint of 
an IMBH. The comparison of the observed mass segregation 
with the simulations clearly demonstrates that the cluster bi- 
nary stars play a key role in defining its dynamical state and 
therefore cannot be ignored. The binary fraction of the cluster 
needs to be better constrained. 
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Fig. 1 . — Map of the entire HST d ata- set used in this work. The position of each camera is labeled. The nominal core radius (r^ = 0.86') and half mass radius 
(rj = 1.81') taken from lHarris et al.l {1996) are showed with dotted and dashed circles respe ctive ly. The full circles define the four annuli in which we divided the 
ACS photometry. For details on the radial division criteria of ACS and WFPC2 see Section [231 
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Fig. 2. — The sharpness (sh) parameter provided as an output by ALLFRAME is used as the source quality diagnostic. By plotting sh as a function of the 
magnitude (left panel) we found that spurious objects, detected around the haloes and along diffraction spikes of saturated stars, have sharpness values that differ 
from the ones representative of bona-fide stars (—0.15 <sh < 0.15). The goodness of the selection in shar pness in identifying false detections is confirmed by the 
fact that the objects with large sh values have also larger photometric errors (right panel; see Section lTTl for details). 
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Fig. 3. — Color magnitude diagram of tlie cluster as obtained from photometric reduction of the WFPC2 (left panel) and ACS (right panel) data sets. The mean 
ridge line as obtained from 2 sigma rejection fitting process along the MS of each catalogue is also showed (open circles). Comparison between the observations 
and models by Baraffe et al. ( 1997) for [M/H]=-1.3 (dashed line) and [ M/H] =-1.0 (full line) is also showed. The shaded regions in the CMDs show the limits 
below which the photometric completeness falls below 50 % (see Section l231 for details). 
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Fig. 4. — Photometric completeness / as a function of I for the two data sets. The ACS and WFPC2FoV have been divided in 4 and 2 area of same completeness 
respectively. The full horizontal line show the limits of 50% of completeness. 
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Fig. 5. — fi/^SMiv OLFs of M 10, divided in different areas (see Section|3j. The measurement as reported in Table|2](full circles in the plot) are shifted by an 
arbitrary amount to make the plot more readable. The theoretical LFs that best fit the data are shown as solid lines. The index of their corresponding power-law 
MF is, from bottom to top, a = 0.7,0.4,0.1,-0.3,-0.6,-0.9. A positive index means that the number of stars decreases with mass. All but the innermost LF 
reach a mass of 0.2Mq . 




Fig. 6. — The LFs of Figure [3] (full circles) ai'e compared with those predicted by our multi-mass Muchie-King model at various radii (reported left of the 
respective LF) inside the cluster (full line), for a GMF with fixed a = —0.7. The predicted radial variation of the LF is fully consistent with the observations. 
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Fig. 7. — Observed radial mass segregation profile (Am(r) measured in Mq) for M 10 (red points with la en'or bars), compared to expectations from numerical 
simulations. The shaded areas in both panels represent the 2a confidence area for our simulations with an IMBH (blue) and without an IMBH (green). The long 
dashed Unes define the la confidence regions. The small inset shows the upper envelope of all snapshots of the simulations including an IMBH in the set. On the 
left panel the simulations contains runs with up to 10% primordial binaries as in Pasauato et al. (2009). On the right panel only the simulations with Miller & 
Scalo IMF, 32k particles and no primordial binaries are use d. Bi naries widen the confidence intervals towards low mass segregation values, thereby preventing a 
firm conclusion as to the presence of an IMBH (see Section lsTl . 




